Conformational Mechanics of Polymer Adsorption Transitions 

at Attractive Substrates 



Monika ModdelQ Michael Bachmann[j] and Wolfhard Jank^l 

Institut fur Theoretische Physik, Universitat Leipzig, Postfach 100 920, D-04 009 Leipzig, 
and Centre for Theoretical Sciences (NTZ), Emil-Fuchs-Strafie 1, D-0 1^105 Leipzig, Germany 

Conformational phases of a semiflexible off-lattice homopolymer model near an attractive sub- 
strate are investigated by means of multicanonical computer simulations. In our polymer-substrate 
model, nonbonded pairs of monomers as well as monomers and the substrate interact via attractive 
van der Waals forces. To characterize conformational phases of this hybrid system, we analyze 
thermal fluctuations of energetic and structural quantities, as well as adequate docking parameters. 
Introducing a solvent parameter related to the strength of the surface attraction, we construct and 
discuss the solubility-temperature phase diagram. Apart from the main phases of adsorbed and 
desorbed conformations, we identify several other phase transitions such as the freezing transition 
between energy-dominated crystalline low-temperature structures and globular entropy-dominated 
conformations. 



I. INTRODUCTION 

The study of the conformational behavior of polymers 
near surfaces is a fascinating field, both from a physical 
and chemical perspective. It provides a rewarding play- 
ground for basic and applied research. With the advent 
of new sophisticated experimental techniques offering an 
enormous potential in polymer and surface research, the 
interest in the hybrid interface of organic and inorganic 
matter has increased. Among such techniques at the 
nanometer scale is, for example, atomic force microscopy 
(AFM), where it is possible to measure the contour length 
and the end-to-end distance of individual polymers [l[ 
or to quantitatively investigate the peptide adhesion on 
semiconductor surfaces Q. Another experimental tool 
with an extraordinary resolution in positioning and ac- 
curacy in force measurements are optical tweezers |3|, H| . 

Applications for adsorption phenomena in polymeric 
solutions can be found in such different fields as lubrica- 
tion, adhesion and surface protection, steric stabilization 
of colloidal particles, as well as in biological processes 
of membrane-polymer interaction. The understanding of 
the latter is particularly important for the reconstruction 
of cell processes. The knowledge of structure formation 
processes near interfaces is also a prerequisite for design- 
ing micro- or nanostructures providing a large variety of 
possible applications in nanotechnology. 

Despite much efforts in the past, the solvent-quality- 
dependent behavior of a dilute polymer solution exposed 
to an adsorbing substrate is not yet fully understood. 
In good solvent, dominating structures are random coils 
since the monomers and the solvent molecules attract 
each other and, consequently, solvent molecules accumu- 
late between monomers and push the monomers apart. 
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Also at high temperatures, random coils are favored as 
they possess a higher conformational entropy than glob- 
ular conformations. Reducing the temperature, the more 
compact low-energy conformations gain thermodynamic 
weight, and the polymer collapses in a cooperative re- 
arrangement of the monomers. At the ^-temperature, 
where the self- avoidance effect and the solvent effect ex- 
actly cancel, the size of a flexible polymer scales like an 
ideal chain, that is, (^g yr ) oc N 2u , with v = 1/2. Globu- 
lar conformations are highly compact but have only little 
internal structure. Hence those conformations are still 
entropy-dominated, and a further transition toward low- 
degenerate crystalline energetic states is expected and 
indeed observed: the freezing transition [5|, |fj]. 

The presence of an attractive surface strongly affects 
the behavior of the polymer in the vicinity of the inter- 
face. The monomer-monomer attraction, being respon- 
sible for the collapse below the #-point, and the surface- 
monomer attraction, resulting in the adsorption, compete 
with each other. This competition gives rise to a variety 
of different conformational phases. The polymer adsorbs 
at the surface, if the temperature is decreased below the 
adsorption transition temperature, but at high tempera- 
tures only a small number of monomers is in contact with 
the surface, even if the polymer was grafted to it @, Hj]- 
This is due to the lower entropy of conformations spread 
out on the surface, compared to the behavior in bulk. 

Numerous detailed studies have been performed to 
elucidate the conformational behavior of homopolymers 
and heteropolymers near substrates. Compared to 
experiments, computer simulations have the advan- 
tage that combinations of parameters can be varied 
at will. Theoretical studies have, for example, been 
performed analytically using scaling theory @, G3, 
mean- field density functional theory PJJ], and series ex- 
pansion [12|, ll3[ and numerically by employing off-lattice 
models such as a bead-spring model of a single polymer 
chain grafted to a weakly attractive surface [H, H3 |. 
multiscale modelling |15| . Monte Carlo studies of lattice 
homopolymers 0, S, Sill, [HI B3, H E3, HH , molecular 
dynamics combined with a stretching of an adsorbed 



homopolymer (22[, or exact enumeration (23[. Also 
adsorption-desorption dynamics were investigated in 
Brownian dynamics simulations of coarse-grained mod- 
els 0. 

In this study, we performed multicanonical Monte 
Carlo computer simulations in order to analyze ther- 
modynamic properties of the adsorption of a semiflexi- 
ble polymer at a flat and unstructured, attractive sub- 
strate. Our main objective is the classification of the 
structural phases accompanying the adsorption process 
and the construction of the complete (pseudo)phase di- 
agram parametrized by the temperature and a suitably 
introduced solvent- quality parameter. The rest of the 
paper is organized as follows. In Sect. HI the hybrid 
polymer-substrate model, the multicanonical simulation 
method, as well as the measured observables are intro- 
duced. The main result, the pseudophase diagram, is pre- 
sented and discussed in detail in Sect. IIII1 Several aspects 
of the phase structure are consolidated by a precise anal- 
ysis of individual observables introduced in Sect. [Ill In 
Sect. II V( our off-lattice results are compared with former 
results obtained in simulations of lattice models. Even- 
tually, in Sect. [VI the paper is concluded by a summary 
of our findings. 

II. THE MODEL 

A. Hybrid Modeling of Polymer-Substrate 
Interaction. 

We employ a coarse-grained off-lattice model for semi- 
flexible homopolymers that has also been generalized for 
studies of heteropolymers [25|, [26[ and helped to under- 
stand protein folding channels from a mesoscopic per- 
spective [27], [28j . In contrast to earlier adsorption stud- 
ies of lattice polymers @, fi i, [H, [HI El EI IH EH , we 
here accept the associated additional computational cost 
of an off-lattice model in order to get rid of undesired 
effects of underlying lattice symmetries. 

As on the lattice, we assume that adjacent monomers 
are connected by rigid covalent bonds. Thus, the distance 
| r^+i — fi\ is fixed and set to unity. Bond and torsional 
angles are free to rotate. The energy function consists of 
three terms, 

E = £?bend + ^LJ + ^sur (1) 

associated with the bending stiffness (i^bend), monomer- 
monomer Lennard- Jones interaction and 
monomer-surface attraction (E sur ). 

A sketch of a coarse-grained polymer segment is de- 
picted in Figure [TJ The position vector of the ith 
monomer, i = 1, . . . , iV, is denoted by r$. A polymer 
with N monomers has N — 1 bonds of length unity be- 
tween neighboring monomers and N — 2 bending angles 
i = 1, . . . , N — 2, defined through 

cos = (r i+ i - n) • (f i+ 2 - n+i) • (2) 
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FIG. 1: A segment of the coarse-grained semiflexible polymer. 
The distance between two adjacent monomers is fixed and 
set to unity. The bending angle at the (i + l)th monomer 
is denoted by and the vector between the ith and jth 
monomer by fij = fj — fi with \ fij \ — rij. 

The Lennard- Jones potential of nonbonded monomers is 
of standard form, 




where = \fj — r^|. The lowest-energy distance of 
the Lennard- Jones potential between two monomers is 
^min = V% ~ 1.12 and is hence slightly larger than the 
unity bond length. The bending energy is defined as 

^ N-2 

E he nd=^ Yl (l-COs(tfi)). (4) 

2=1 

The angle di is in the interval [0, tt) and the bending en- 
ergy can be considered as a penalty for successive bonds 
deviating from a straight arrangement. 

The attractive surface potential is obtained by in- 
tegrating over a smooth half-space, where every space 
element interacts with a single monomer by the usual 
Lennard- Jones 12-6 expression. One obtains the poten- 
tial [H,^: 




where Z{ is the distance of the ith monomer from the 
surface. The parameter e s defines the surface attraction 
strength, and as such it weighs the energy scales of intrin- 
sic monomer-monomer attraction and monomer-surface 
attraction. 

A distance z = Lb ox away from the attractive surface, 
we place a steric wall that is necessary to prevent the 
molecule from escaping. Because the exact form of the 
density of states depends on the box height Lb ox , all mea- 
sured quantities also depend on the choice of Lb ox - As 
soon as the box size exceeds the polymer size, however, 
the influence on the observables, apart from the substrate 
distance of the center-of-mass z cm = ^2iLi %i/N of the 
polymer, is sufficiently small. We chose Lb ox = 20 (40) 
for the polymer with N = 13 (20) monomers. 



B. Multicanonical Method. 

The canonical partition function of our hybrid system 
at temperature T is given in natural units by 

Z= / dEg(E)e- E / T , (6) 

where g(E) = e s ^ is the density of states that con- 
nects (microcanonical) entropy S and energy. Therefore, 
all information regarding the phase behavior of the sys- 
tem - typically governed by the competition between en- 
tropy and energy - is encoded in g(E). Consequently, 
for a detailed global analysis of the phase behavior, a 
precisely estimated density of states is extremely helpful. 
Unfortunately, the density of states covers many orders 
of magnitude in the phase transition regimes, so that 
its estimation requires the application of sophisticated 
generalized-ensemble Monte- Carlo methods. 

For our analyses, we have performed multicanonical 
simulations [3l|, as multicanonical sampling allows the 
estimation of g(E), in principle, within a single simula- 
tion. The idea of the multicanonical method is to increase 
the sampling rate of conformations being little favored 
in the free-energy landscape and, finally, to perform a 
random walk in energy space. This is achieved in the 
simplest way by setting T = oo and introducing suitable 
multicanonical weights W muca (E) ~ g~ 1 (E) in order to 
sample conformations X according to a transition prob- 
ability 

u(X -> X') = min[e s ( B ( x »- s ( £ ( x '>>, 1], (7) 

where S(E(X)) = - In W muca (E(X)) = lng(E(X)). 

The implementation of multicanonical sampling is not 
straightforward as the multicanonical weights W muca (E) 
are obviously unknown a priori. Therefore, starting with 
W^ uca (^) = const., the weights have to be deter minded 
by an iterative procedure until the multicanonical his- 
togram is almost "flat", that is, if the estimate for the 
density of states after the nth run, g^ n \E), satisfies 

0< n >(25)W^(25) « const. (8) 

in the desired range of energies. An efficient, error- 
weighted estimation method for the multicanonical 
weights is described in detail in refs [32[ and [33( | . 

Eventually, if eq [8] is reasonably satisfied, the multi- 
canonical weights ^muca(^) = [9^ n \E)]~ l are then used 
in a final long production run, where all quantities of in- 
terest are measured and stored in a time-series file. The 
canonical expectation value of any quantity O at temper- 
ature T is then obtained from the multicanonical time 
series of length M by reweighting, 

(0) = Sgi Q(Xt)e- E ^» T W-UE(Xt)) ( ) 

E M ie - B (XO/T^- u l ca( £ (Xt)) ' 

where t is the multicanonical Monte Carlo "time" step 
(or sweep). 
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FIG. 2: Phase diagram of the 20mer. The colored stripes 
separate the individual conformational phases (see text and 
Figure [3} . 



C. Suitable Energetic and Structural Quantities for 
Phase Characterization. 

To obtain as much information as possible about the 
canonical equilibrium behavior, we define the following 
suitable quantities O. Next to the canonical expecta- 
tion values (O) , we also determine the fluctuations about 
these averages, as represented by the temperature deriva- 
tive ({OE) — (O) (E)) /T 2 . We use generic units, in 
which ks = 1- 

Apart from energetic fluctuations such as the specific 
heat, cy = d (E/N) dT, several structural quantities are 
of interest. The radius of gyration is a measure for 
the extension of the polymer and defined by R 2 yY = 

ElAifi - Rcm) 2 )/N = E£iEjLi((ri - r 3 ?)/2N* 

with R cm = ^f =1 Ti/N being the center-of-mass of the 
polymer. Since the substrate introduces a structural 
anisotropy into the system, it is not only useful to inves- 
tigate the overall compactness of the polymer expressed 
by (R gyr ), but also to study the expected different behav- 
ior of its components parallel and perpendicular to the 
surface: R\ = £ti Ejli<(*i " *jf + {Vi ~ Vi?)/™ 2 
and R\ = Eti EjLi((*i " Zjf)/2N 2 such that R 2 gyl = 

R 2 + R\ . 

Clear evidence that the polymer is, on average, freely 
moving in the box or very close to the surface can be 
provided by the distance of the center-of-mass, z cm , of 
the polymer to the surface. 

Another useful quantity is the mean number of 
monomers docked to the surface. A single-layer structure 
is formed if all monomers are attached at the substrate; 
if none is attached, the polymer is quasifree (desorbed) 
in solvent. The surface potential is a continuous poten- 
tial, and in order to distinguish monomers docked to the 
substrate from those not being docked, it is reasonable 
to introduce a cutoff. We define a monomer i as being 
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pseudophase typical configuration 
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FIG. 3: Representative examples of conformations for the 
20mer in the different regions of the T-e s pseudophase dia- 
gram. DE, DG, and DC represent bulk "phases", where the 
polymer in preferably desorbed. In regions AE1, AE2, AC1, 
AG, AC2a, and AC2b, conformations are favorably adsorbed. 



"docked" if zi < z c = 1.2. The corresponding measured 
quantity is the average ratio (n s ) of monomers docked to 
the surface and the total number of monomers. This can 
be expressed as n s = N s /N with N s = ©(^c — Zi), 

where Q(z) is the Heaviside step function. 

Similary, the mean number of intrinsic contacts 
is introduced as a measure of the global compact- 
ness of the polymer: n rn = N m /N with N m = 
E^ 7 IT 2 E^L^+2 e ( e c-eLJ(rij)), where e c = -0.2 and 
e LJ (^)=4(r-. 12 -r- 6 ). 




FIG. 4: (a) Specific heat as a function of temperature T and 
surface interaction strength e s for the 20mer. Lines represent 
the simulation data for fixed values of e s , the color code is 
interpolated, (b) Specific heat curves for different values of 



sweeps and was performed with at least 2 different ini- 
tializations of the random number generator in order to 
avoid systematic deviations. 

For a convenient overview, we display our final pseu- 
dophase diagram of the 20mer already here in Fig. [2j 
Conformations representative of the different phases are 
shown in Fig. OH The details will be discussed in the fol- 
lowing. It should be stressed that all phases and transi- 
tions occuring in our analysis are not phases in the strict 
thermodynamic sense, since we are dealing with finite 
chain lengths. However, even for the rather short chains 
considered here, a reasonable picture of the polymer ad- 
sorption behavior at the surface is obtained, and most of 
the phases are believed to still exist for longer chains. 



III. RESULTS AND DISCUSSION 

Our major objective is the construction of the pseu- 
dophase diagram of conformational phases, based on 
our results obtained for energetic and structural fluctu- 
ations. To this end, multicanonical simulations [3l| for 
51 different surface attraction strengths e s , in the range 
e s = 0, ...,5, were performed and reweighted to tem- 
peratures T G (0,5]. Each simulation consisted of 10 8 



A. Energetic Fluctuations 

In Fig. Ufa), the specific heat of the 20mer is shown 
as a function of T and e s . Although for both investi- 
gated polymer lengths the energy varies smoothly with 
changing T and e s (with a global minimum at minimal 
temperature and maximal surface attraction), two tran- 
sitions can be identified as ridges in the profile of the 
specific heat: The first one is the adsorption transition 
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separating desorbed and adsorbed conformations. The 
other transition is sort of a freezing transition at low 
temperatures. Near T = 0.25, the specific heat exhibits a 
pronounced peak and decreases rapidly at lower temper- 
atures, independently of the surface attraction strength. 
The low-temperature peak in specific heat and the crys- 
talline shape of structures found in this regime signalize 
this freezing transition. Although the freezing tempera- 
ture seems to be rather constant, the type of crystalline 
structure adopted by the polymer depends strongly on 
e s . For the identification of the polymer shapes, we take 
in the next subsection a closer look at the conformational 
quantities. 




B. Comparative Discussion of Structural 
Fluctuations 

The radius of gyration provides an excellent measure 
of the globular compactness of polymer conformations. 
Figures [5] and [6] reveal that the most compact conforma- 
tions dominate at low temperatures and for small values 
of surface attraction strengths e s . The freezing transition 
can be found at the same temperature, as it has already 
been identified from the specific heat. The adsorption 
transition typically affects only segments of the polymer 
and is thus not prominently signaled by the radius of 
gyration. 

This is different for its components parallel (Figures [3 
and[8j) and perpendicular (Figures l9l and [TO]) to the sur- 
face, respectively. For example, for e s > 3.4, (R±) van- 
ishes at low temperatures for the 20mer, whereas (R\\) 
attains low values at lower e s . The vanishing of (R±) cor- 
responds to conformations where the polymer is spread 
out flat on the surface without any extension into the 
third dimension. The associated pseudophases are called 
adsorbed compact (AC1) and adsorbed expanded (AE1) 
phase. The '1' is appended in order to distinguish these 
phases from topological three-dimensional phases, such 
as, for example, the AC2 subphases. The phases AC1 
and AE1 are separated by the freezing transition such 
that polymer structures in AC1 are maximally compact 
at lower temperatures, whereas AE1 conformations are 
less compact and more flexible but still lie rather planar 
at the surface. To verify that conformations in AC1 
are indeed maximally compact single layers, we argue as 
follows. The most compact shape in the two-dimensional 
(2D) continuous space is the circular disk. Thus, one can 
calculate (R\\) for a disk and compare it with the sim- 
ulated value. Assuming N monomers to be distributed 
evenly in the disk, N « 7rr 2 , where r is the radius of 
the disk in units of the mean distance of neighboring 
monomers. The radius of gyration in the same units is 
thus given by 

Since we have two types of mean distances between 




FIG. 5: (a) (Rgyr) of the 20mer as a function of T and e s . (b) 
(Rgyr) as a function of the temperature for various values of 
e s . 

monomers in compact conformations dependent on 
whether they are adjacent along the chain or not, 
w e expe ct for disk-like conformations an the s urface : 
^20/2^ « 1.784 < (#11,20) < 2.003 « r min ^20/2^. 
The simulated value is (#11,20) ~ 1.81, which nicely fits 
the estimate. The equivalent estimate for the 13mer also 
confirms discoidal conformations in AC1. 

The argument is similar for sphere-like three- 
dimensional (3D) compact conformations with TV = 
47rr 3 /3. Corresponding conformations are found as free 
desorbed compact chains (DC), as well as adsorbed com- 
pact polymer conformations (AC2a) for weak surface at- 
traction. In this case, the radius of gyration is given by 

2 1 f ^3 / / 2 3 2 3 {3N\ 2/S 

(ii) 

The estimates (3/5) 1 / 2 (3 x 13/47T) 1 / 3 « 1.130 < 
(# g yr,i3> < 1-268 * r min (3/5) 1 / 2 (3 x 13/47T) 1 / 3 , and 
1.684 < (i? gy r,2o) < 1.464 slightly overestimate the simu- 
lated values (iigyr,i3> = 1-023 and (# gyr ,20> = 1-242. The 
deviations can be explained by the fact that the mass of 
the polymer is not uniformly distributed in the sphere as 
it is assumed in the calculation. For a compact packing 
of discrete monomer positions, it is more realistic that 
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FIG. 6: (a) d (R gyr ) /dT of the 20mer. (b) d {R gyr ) /dT, 
parametrized by e s . 



the outer thin shell of the sphere does not contain any 
monomers. Performing the integration not from r' = 
to r' = r, but only to r' = r — e, reduces the estimated 
radius of gyration significantly already for small e due to 
the increased weight of the outer shells in higher dimen- 
sions. Taking this effect into account, the thus obtained 
values of (R gy r) seem to be even more reasonable. 

The most pronounced transition is the strong layering 
transition at e s « 3.4 for N = 20 that separates regions 
of planar conformations (AC1, AE1) in the T-e s plane 
from the region of stable double-layer structures (AC2b) 
and adsorbed globules (AG), below and above the freez- 
ing transition, respectively. For high surface attraction 
strengths e s , it is energetically favorable to form as many 
surface contacts as possible. In the layering-transition 
region, a higher number of monomer-monomer contacts 
causes the double-layer structures to have just the same 
energy as single-layer structures. For lower e s values, 
the double-layer structures possess the lowest energies. 
Hence this transition is a sharp energetical transition. 

Although for the considered short chains no higher- 
layer structures are observed, the components (Rn j_) in- 
dicate some activity for low surface attraction strengths. 
For N = 20, e s « 1.4 is the lowest attraction strength, 
where still stable double-layer conformations are found 
below the freezing transition. What follows for lower e s 
values after a seemingly continuous transition is a low- 
temperature subphase of surface attached compact con- 
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T 

FIG. 7: (a) (R\\) of the 20mer. (b) (R\\) for selected values 
of e s . 




0.5 1 1.5 2 2.5 



FIG. 8: (a) d(R {l ) /dT for the 20mer. (b) d(R {l ) /dT for 
different surface attraction strengths e s . 
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formations, called AC2a. AC2a conformations occur if 
the monomer-surface attraction is not strong enough to 
induce a layering in compact attached structures. The 
characterization of structures in this subphase requires 
some care, as system-size effects are dominant. Although 
the surface attraction is sufficiently strong to enable 
polymer-substrate contacts, compact desorbed polymer 
conformations below the ^-transition are not expected to 
change much. Thus, layering effects do not occur. We 
found two distinct classes of structures in this region: (1) 
completely undistorted compact conformations touching 
the surface and (2) semispherically shaped structures 
docked at the surface. For the rather short chains in our 
study, the occurrence of these shapes strongly depends 
on the precise number of monomers. Not surprisingly, 
we find differences for the 13mer and the 20mer. 

For N = 13, both components of the radius of gyra- 
tion, d (R±) /dT and d (Ru ) /dT, indicate a transition at 
e s « 0.45. Low-energy configurations reveal that this is 
a wetting transition between undistorted compact con- 
formations for smaller e s and docked conformations for 
largerr surface attraction strengths. An analogous tran- 
sition for N = 20 was not found. In this case, the AC2a 
pseudophase seems to consist of a mixture of both types 
without any transition between them. This is also con- 
firmed by analyzing the low-energy conformations in this 
regime. The higher the e s value, the larger the average 
number of monomers being docked at the surface, but a 
clear cut from the compact adsorbed conformations does 
not exist. This difference in the wetting transition for 
N = 13 and N = 20 might be due to the fact that the 
most compact conformation for N = 13 is an almost per- 
fect icosahedron [34| . "Almost" because the length scales 
of covalent bonds and intermonomeric Lennard- Jones in- 
teraction differ slightly. 

We also searched for low-energy states with a modified 
LJ energy minimum distance shifted to unity and indeed 
found perfectly icosahedral morphologies. This addition- 
ally stabilizes the polymer conformation and is already 
known from studies of atomic clusters. The smallest 
Mackay icosahedron [35[ with characteristic 5-fold sym- 
metry is formed by 13 atoms. Larger perfect icosahe- 
dra also require a "magic" number (55, 147, 309,. . . ) 
of atoms. This holds also true for crystals of elastic 
polymers |34j. Thus, it might be worthwhile to also study 
the wetting transition for other chain lengths in order to 
be able to predict a trend for longer chains, which is 
not possible only knowing the behavior for the two chain 
lengths investigated in our study. The parameters of the 
low-temperatures pseudophase transitions for the 13mer 
and the 20mer are summarized in Table H 



Raising the temperature above the freezing tempera- 
ture, polymers from adsorbed and rather compact con- 
formations that look like glocular, unstructured drops on 
the surface. This pseudophase is called a surface- attached 
globule (AG) phase. It has been first conjectured from 



TABLE I: Surface Attraction Strength e s for all Low- 
temperature Transitions for N = 13 and N = 20 
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e s 


« 0.2 


transition AC «-»■ AC2a 


e s 


« 0.5 






transition AC2a AC2b 


e s 
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layering transition AC2b <-> AC1 


€s 


« 2.8 


€s 


& 3.4 




FIG. 9: (a) (R ± ) of the 20mer. (b) (R ± ) for selected values 

of € 8 . 



short exact enumeration studies of 2D polymers in poor 
solven t |23| but was also found in lattice-polymer stud- 
ies @, Il6|. At even higher temperatures, two scenar- 
ios can be distinguished in dependence of the relative 
strengths of monomer-monomer and monomer-substrate 
interactions. In the first case, the polymer first desorbs 
from the surface [from AG to the desorbed globular (DG) 
bulk phase] and disentangles at even higher temperatures 
[from DG to the desorbed expanded bulk phase (DE)]. 
In the latter case, the polymer expands while it is still 
on the surface (from AG to AE2) and desorbs at higher 
temperatures (from AE2 to DE). Due to the higher rel- 
ative number of monomer-monomer contacts in compact 
bulk conformations of longer chains, the ^-temperature 
increases with N . The same holds true for the surface 
attraction strength, e s , associated with the layering tran- 
sition. 
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FIG. 10: (a) d(R±) /dT of the 20mer. 
different surface attraction strengths e s . 



(b) d(i^_L> /dT for 
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FIG. 11: (a) (2cm) °f the 20mer. (b) (z cm ) for various values 
of the parameter e s . 



FIG. 12: Mean number of surface contacts (n s ) vs. surface 
attraction strength e s for the 13mer and the 20mer at T = 
0.001. 



It is clear that the structural behavior of the studied 
small chains is affected by finite-size effects, in particular 
in the compact pseudophases. As long as surface effects 
are as influential as volume effects, the shapes of com- 
pact adsorbed (but also of compact desorbed 34]) con- 
formations differ noticeably for polymers with different 
but small lengths, and a precise classification is difficult. 
However, for longer chains, filmlike (AC1) and semispher- 
ical conformations (AC2), as well as surface-attached 
droplets (AG), will dominate the respective phases. Cur- 
rently, the simulation of longer chains, aiming at the iden- 
tification of all conformational phases and the quantita- 
tive analysis in the thermodynamic limit, is too chal- 
lenging. Thus, a more detailed classification within the 
compact phases is left for future work. 



C. Adsorption Parameters 

The adsorption transition can be discussed best when 
looking at the distance of the center-of-mass of the poly- 
mer to the surface (Figure ITT]) and the mean number of 
surface contacts (Figures [12] and [13]) . As can be seen 
in Figure [HI for large temperatures and small values of 
e s , the polymer can move freely within the simulation 
box and the influence of the substrate is purely steric. 
Thus, the average center-of-mass distance (z cm ) of the 
polymer above the surface is just half the height of the 
simulation box. On the other hand, for high e s values 
and low temperatures, the polymer favors surface con- 
tacts and the average center-of-mass distance converges 
to (z cm ) « 0.858, corresponding to the minimum-energy 
distance of the surface attraction potential for single- 
layer structures, and slightly larger valuesr for double- 
layer and globular structures. 

One clearly identifies a quite sharp adsorption tran- 
sition that divides the projection of (z cm ) in Figure 
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FIG. 13: (a) (n s ) of the 20mer. (b) (n s ) for different values 
of e s . 



[TTli into an adsorbed (bright /green) regime and a des- 
orbed (dark/blue) regime. This transition appears as a 
straight line in the phase diagram and is parametrized 
by e s oc T. Intuitively, this makes sense since at higher 
T the stronger Brownian fluctuations are more likely to 
overcome the surface attraction. 

Consistently with our above discussion, (z cm ) and 
d(^cm) /dT also reveal the low-temperature transitions 
between the adsorbed phases AC2a, AC2b, and AC1. 
for the detailed discussion of these adsorbed phases, we 
concentrate ourselves on the mean number of surface con- 
tacts (n s ) (Figs. 

Unlike in simple-cubic lattice studies, where one finds 
(n s ) ~ 1 for an /-layer structure [7[, we find for double 
layer structures (n s ) > 1/2. The reason is that most 
compact multilayer structures are cuboids on the lattice, 
whereas in our off-lattice study, "layered" conformations 
correspond to semispherical shapes, where, for optimiza- 
tion of the surface of the compact shape, the surface layer 
contains more monomers than the upper layer. Since this 
only regards the outer part of the layers, the difference 
is more pronounced the shorter the chain is. 

In Figure \12\ (n s ) is shown as a function of e s for small 
temperatures. (n s ) is a useful quantity to identify lay- 
ering effects. Starting at high e s , for both chain lengths 
first (n 3 ) « 1 until at the layering transition, (n s ) jumps 



to (n s ) « 0.69 for N = 13 and to (n s ) « 0.65 for N = 20. 

Further jumps corresponding to further layering transi- 
tions are not observed. Instead, what follows is a plateau 
regime where the relative amount of monomers that cover 
the surface is rather constant. When the double-layer 
structure gets unstable at lower e s , (n s ) starts to de- 
crease again. The conformations in AC and AC2a thus 
do not exhibit a pronounced number of surface contacts, 
and (n s ) varies with e s . Near e s ~ 0.2, where the polymer 
desorbs, (n s ) converges rapidly to (n s ) = as e s — > 0. To 
conclude, the single- to double-layer "layering transition" 
is a topological transition from 2D to 3D polymer confor- 
mations adsorbed at the substrate. The solvent-exposed 
part of the adsorbed compact polymer structure, which 
is not in direct contact with the substrate, reduces under 
poor solvent conditions the contact surface to the solvent. 
Because of the larger number of degrees of freedom for 
the off-lattice polymer, layered structures are not favored 
in this case. Thus, higher-order layering transitions are 
not identified in our analyses (which in part is also due to 
the short lengths of the chains) but are also not expected 
in pronounced form. 

The observable left to discuss is the mean number of 
intrinsic contacts. It behaves very much like the radius 
of gyration, such that the projection of (n m ) onto the 
T-e s -plane is divided into a compact regime comprising 
AC, AG, AC2a, AC2b, DC, and DG and a regime of less 
compact conformations. This nicely confirms the results 
already obtained. 



D. The Pseudophase Diagram 

To summarize all the informations gained from the dif- 
ferent observables, we construct the approximate bound- 
aries of different regimes in the T-e s plane. The pseu- 
dophase diagram was already displayed in Figure [2] where 
the different pseudophases are denoted by the abbrevia- 
tions explained in the previous subsection. 

The pseudophases found are (for selected representa- 
tive conformations see Figure [3]): 

• DE (desorbed expanded): Random-coil phase of 
the quasifree desorbed polymer. 

• DG (desorbed globular): Globular phase of the 
desorbed chain. 

• DC (desorbed compact): Maximally compact, 
spherically shaped crystalline structures dominate 
this desorption phase below the freezing-transition 
temperatue. 

• AE1 (adsorbed expanded single layer): Adsorbed 
phase of expanded, rather planar but little compact 
random-coil conformations. 

• AE2 (adsorbed expanded 3D conformations): Ad- 
sorbed, unstructured random-coil-like expanded 
conformations with typically more than half of the 
monomers in contact with the attractive substrate 
are favored in this pseudophase. 
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• AC1 (adsorbed compact single layer): Phase of ad- 
sorbed circularly compact filmlike conformations. 

• AG (adsorbed globular 3D conformations): Repre- 
sentative conformations are surface-attached glob- 
ular conformations and look like drops on the sur- 
face. 

• AC2a (adsorbed compact 3D conformations): 
Compact, semispherically shaped crystalline con- 
formations dominate in this subphase. 

• AC2b (adsorbed compact double layers): Sub- 
phase of adsorbed, compact double-layer conforma- 
tions. The occupation of the surface layer is slightly 
larger than that of the other layer. 

AC2a and AC2b are subphases in the regime of the phase 
diagram, where adsorbed compact and topologically 
three-dimensional conformations are dominant. Since 
pronounced layering transitions as observed in lattice- 
polymer studies are not expected here, the discrimina- 
tions of AC2a and AC2b is likely to be irrelevant in the 
thermodynamic limit. However, AC2a,b differ qualita- 
tively from the phase AC1 of topologically 2D polymer 
films and we thus expect that the transition between film- 
like (AC1) and semispherical conformations (AC2) is of 
thermodynamic relevance. 

For N = 13, where the maximally compact confor- 
mation is more stable due to the high symmetry of the 
icosahedral structure, we found an additional subphase: 

• AC (adsorbed icosahedral compact conforma- 
tions): Like DC, but polymers typically are in 
touch with the surface. As a clear individual sub- 
phase only observed for N = 13. 

The transition lines in the pseudophase diagram (Fig- 
ure [2]) represent the best compromise of all quantities 
analyzed separately in our study. Only in the thermo- 
dynamic limit of infinitely long chains are most of the 
identified pseudophase transitions expected to occur at 
sharp values of the parameters e s and T for all observ- 
ables. For finitely long chains, the transition lines still 
vary with chain length N and are not well defined due to 
the broad peaks that are slightly different for different ob- 
servables (see Figure [T4|). Taking that into account, the 
pseudophase diagram gives a good qualitative overview 
of the behavior of polymers near attractive substrates 
in dependence of environmental parameters such as sol- 
vent quality and temperature. The locations of the phase 
boundaries should be considered as rough guidelines. 



IV. COMPARISON WITH LATTICE RESULTS 

Finally, we would like to compare our results discussed 
here with those obtained from a similar model on a 
simple-cubic (sc) lattice pH [l8|. 

The lattice polymer is modeled as a nongrafted inter- 
acting self-avoiding walk confined between two infinitely 
extended parallel planar walls. One wall is short-range 
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FIG. 14: (a) Temperature dependence of several observables 
for e s = 2 and N = 20. (b) The derivative with respect to T 
of the same quantities. 



attractive, and the other is purely sterical and prevents 
the polymer from escaping. The energy of the system is 
given by 



e L n L 

m mi 



(12) 



where n\ is the number of nearest-neighbor monomer- 
substrate contacts, is the number of nearest-neighbor 
but nonadjacent monomer- monomer contacts, and 
and are the respective contact energy scales. In 
refs [16] and [18] the contact density g ns m m was di- 
rectly sampled by means of the contact-density chain- 
growth method, which is an extension of the multicanon- 
ical chain-growth method [36|, [37| • The pseudophase 
diagram, parameterized by temperature and monomer- 
monomer interaction strengh was discussed mainly using 
the specific-heat profile. For a review see ref [38]. The 
surface-monomer attraction strength was fixed. With the 
contact density, the specific-heat profile can be calcu- 
lated for fixed monomer- monomer interaction = 1, 
while varying the surface attraction parameter as it 
was done in the present off-lattice study. The resulting 
pseudophase diagram is depicted for a 179mer in Fig- 
ure [T5J where the parametrization chosen in ref [l6[ was 
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FIG. 15: Pseudophase diagram of a lattice polymer with 
179 monomers[16], parametrized by the surface attraction 
strength, e^, and temperature, T. The color encodes the 
specific-heat profile; the darker the color, the larger its value. 

rescaled in order to allow for a more direct comparison 
with the results of our off-lattice study. Denoting the 
energy and temperature from ref [16] with E' and T", 
respectively, the rescaling works as follows: 

(13) 



Here, s = e^/e^ is the ratio of energy scales of intrinsic 
and surface contacts as introduced in ref fl6| . It can be 
interpreted as a reciprocal solvent-quality parameter. 

Certain similarities between the off-lattice (Figure [2]) 
and the lattice pseudophase diagram (Figure [T5|) are ob- 
vious. For instance, the adsorption transition line is 
parametrized in both models by e s ex T. Different, how- 
ever, is not only the slope, which depends on the system's 
geometry and energy scales, but also for the off-lattice 
model the extrapolation of the transition line seems to 
go through the origin e s = and T = 0, whereas there is 
an offset observed in the lattice-system analsis such that 
the extrapolated transition line roughly crosses = 0.4 
and T = 0. We speculate that this might be due to 
the intrinsic cuboidal structure of the polymer confor- 
mation on the sc lattice that possess planar surfaces at 
low temperatures even in bulk. Unlike for off-lattice 
models, where a compact polymer attains a spherical 
shape, such a cuboidal conformation is likely to dock at a 
substrate without substantial conformational rearrange- 
ments. Here lies an important difference between lattice 
and off-lattice models. The off-lattice model provides, 
for sufficiently small surface attraction strengths, a com- 
petition between most compact spherical conformations 
that do not possess planar regions on the polymer surface 
and less compact conformations with planar regions that 
allow for more surface contacts but reduce the number of 



intrinsic contacts. 

This also explains why a transition like the one ob- 
served for N — 13 between AC and AC2a, the wetting 
transition, is more difficult to observe in adsorption stud- 
ies on regular lattices. On the other hand, AC2 confor- 
mations at low T and for e s between the adsorption and 
the single-double layering transitions can be observed in 
both models. Similarly in both models, there exists the 
AG pseudophase of surface-attached globules. 

Whereas for the off-lattice model, apart from the wet- 
ting transition, there is only the transition from AC2a 
(semispherically shaped) to AC2b (double-layer struc- 
tures), on the lattice AC2 comprises of a zoo of sub- 
phase transitions. These are higher-order layering tran- 
sitions. Decreasing the surface attraction at low tem- 
perature, layer after layer is added until the number of 
layers is the same as in the most compact conformation. 
A lattice polymer has no other choice than forming lay- 
ers in this regime. The layering transition from AC1 
to AC2 is very sharp in both models. Also the shape 
of the transition region from topologically 2D adsorbed 
to 3D adsorbed conformations looks very similar. Inter- 
estingly, the e^/e^-ratio predicted for this transition in 
ref [7[ agrees quite well with that observed in our off- 
lattice study. For low-energy conformations, it is ar- 
gued that I 3 / 2 = (1 — e\ /e^jA/ -1 / 2 on the square lattice. 
With I = 1.5 and N = 179 this gives e^/e^ = 0.914 
for the single- to double-layer transition, which is con- 
firmed by Fig. [151 We re-expressed this argument for a 
triangular lattice, which describes the low-temperature 
conformations of our off-lattice model better, and ob- 
tain Z 3 / 2 = 2(3 — e s /e m )Af 1//2 /3 for low-energy configura- 
tions. This yields the larger ratio e s /e m = 2.235 (2.384) 
for N = 13 (20), which is due to the higher coordination 
number of this geometry. It is in good agreement with 
our simulation results for this transition. The higher co- 
ordination number also causes the different slopes of the 
respective adsorption transitions. 

To summarize, we conclude that, in particular, the 
high-temperature pseudophases, DE, DC/DG, AG, and 
AE, nicely correspond to each other in both models. 
Noticeable qualitative deviations occur, as expected, in 
those regions of the pseudophase diagram where compact 
conformations are dominant. 



V. SUMMARY 

In this paper, we have constructed the pseudophase 
diagram of thermodynamic conformational phases of a 
single semiflexible homopolymer near an attractive sub- 
strate in dependence of the external parameters surface 
attraction strength and temperature. 

For two polymer chains with N = 13 and N = 20 
monomers, respectively, the canonical expectation val- 
ues of several energetical and structural quantities and 
their thermal fluctuations were measured in multicanon- 
ical computer simulations over a broad range of surface 
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attraction strengths and temperatures. Conformational 
phases and phase boundaries and their location in the 
pseudophase diagram were identified in precise analyses 
of structural fluctuations and suitable adsorption param- 
eters. 

Although the computational expense to accurately ex- 
plore such a broad parameter range restricted us to inves- 
tigate rather short chains, we identified the following con- 
formational pseudophases and pseudophase transitions: 

• Crystalline structures in the regimes of compact 
phases. We identified maximally compact des- 
orbed conformations in bulk (DC) or adsorbed at 
the substrate (AC), semispherical compact confor- 
mations (AC2a) that are distorted by the surface 
but not layered, double-layer conformation (AC2b), 
and single-layer conformations (AC1). 

• Adsorbed conformations in the globular and ex- 
panded (random-coil) phases. Here, three con- 
formational pseudophases were distinguished: un- 
structured 3D surface-attached globular conforma- 
tions (AG), expanded dissolved but planar ad- 
sorbed conformations (AE1), and 3D expanded 
random-coil-like adsorbed conformations (AE2). 

• Desorbed conformations. Compact conformations 
(DC) are separated by the freezing transition from 
globular conformations (DG). At even higher tem- 
peratures above the ^-transition, random-coil con- 
formations are found (DE). 

The sharpest pseudophase transition identified is the 
layering transition between single- and double-layer- 
structures. Higher-layer conformations were not found 



for these short chains. Unlike in recent studies on a 
simple-cubic lattice, where for weak surface attraction 
and positive self-attraction, layering transitions were ob- 
served until a maximally compact cubic structure is 
reached, off-lattice polymers favor maximally compact 
spherical conformation. Thus, we find the expected dif- 
ferences in the behavior of off-lattice and lattice polymers 
in phases, where compact adsorbed conformations dom- 
inate. For the majority of pseudophases, in particular 
those that are assumed to be relevant in the thermody- 
namic limit, we find, however, a nice qualitative coinci- 
dence. This similarity demonstrates the ability of such 
simple coarse-grained models to capture the general ad- 
sorption behavior of polymers near attractive surfaces. 
The increasing experimental and technological capabili- 
ties should allow not only for the experimental verifica- 
tion of the described thermodynamic phases but also for 
a detection of the pseudophases of finite polymers. Since 
polymers are naturally of finite length, this problem is 
one of real physical interest. 
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